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Abstract 

In this article, we provide two different ways of producing recursive formu- 
las for the Rayleigh functions of eigenvalues of the one-dimensional Laplacian 
with non-local boundary conditions, which commutes with an integral opera- 
tor having a harmonic kernel. 
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1. Introduction 

The purpose of this article is to prove recursion formulas for power sums 
for the eigenvalues of a problem which has arisen in recent work by one of us 
Gl. We concentrate on the one-dimensional case for the non-local boundary 
value problem (BVP) described in [1]. Our recursion formulas emulate those 
developed by various authors for Rayleigh functions, or power sums, involving 
roots of various transcendental equations. It was Euler who first found the first 
few closed expressions for what later came to be known as the Rayleigh function 
Q (see also |3, Sec. 15.5], fl): 
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where j Vi „ denotes the n-th positive root of z~ v J v (z), and J v (z) is the Bessel 
function of the first kind of order v. Euler's method was further developed by 
Lord Rayleigh jstl and Carlitz J(l . Both Euler and Rayleigh analyzed eigenval- 
ues of oscillations of physical systems (a hanging chain for Euler and a circular 
membrane for Rayleigh), which aroused their interest in computing zeros of 
the Bessel functions. By exploiting a differential equation of Riccati-type satis- 
fied by the function z~ v J v {z), Kishore developed recursion formulas for 
(Jieiy), starting with the known expression, due to Euler and Rayleigh 
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In his famous book [ 10] , Rayleigh was further led, in the context of treating 
the transverse vibrations of a clamped beam, to finding summation formulas 
for the reciprocal 4th and 8th powers of the positive roots of the equation 



cos x cosh x± 1 = 0. 



(2) 



The early history of the techniques of proving these power sum rules can 



be found in Watson's book [3] as well as till [12J, [13J . The more recent articles 



1 14 , 15 , 3 offer modern views, survey recent results, and apply the techniques 
to various transcendental functions. 

Properly speaking, the technique of resolution of many of these problems 
goes back to Euler and his famous resolution of the "Basel" problem, named af- 
ter the native Swiss city of Euler and the Bernoulli brothers. Euler successfully 
solved the problem first posed by Pietro Mengoli in 1644 (T2L [HI and found a 
closed form for the expression L^ =1 It is now folklore that the sum is 
Heuristically, Euler's argument of 1740 1 17] (see also lL2lll3|]) amounted to writ- 
ing in two different ways: as a Maclaurin series and as the infinite product 
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since the roots of the transcendental equation sin x/x = are given by x = ± nn, 
for n = 1,2, .. . Expanding the product, and equating the coefficients of x 2 gives 
the above formula. For rigorous justifications of these formulas one should 
consult [18, Chap. 1]. Euler's technique is exactly what Rayleigh employed in 
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the case of equation fl2). Many nice examples illustrating this technique appear 
in the excellent paper of Speigel Hil where generalizations of Newton's known 
formulas for the symmetric sums of the roots of a polynomial can be found (see 
also the comments in II al l. 

Radoux 1 20], Liron 121, 22, 23], and more recently Gupta-Muldoon 1 14] and 



Ismail-Muldoon fill employed similar techniques to generate various recur- 
sion formulas in the same spirit. In the case of Radoux and Liron, one finds 
explicit and recursive formulas for sums of even powers of reciprocals for the 
roots of the equation tanx = x, and cotx = x. To illustrate the case of the equa- 
tion, tanx = x, with x\, X2, . . . denoting the strictly positive roots of the equation, 
they derived the sums of even powers of x^'s, i.e., x k , p = 1,2,.. .. For ex- 
ample, the cases p = 1,2 lead to 



OO 1 1 

y — = — 



OO 1 1 

jr_L=_L 

fc=l x k 350 

All of these are manifestations of convolution formulas relating the trace of 
the compact operator defined by the Green's function, and power sums of the 
eigenvalues as detailed i n 12 411 and the classical book of Mikhlin [25] . The recent 
survey paper of Grieser [26] offers a view that relates these formulas to what is 
known for matrices. 

Radoux |20] attributes the method of finding sums of reciprocals of powers 



of eigenvalues of certain operators to Serge Nicaisse, but as detailed in |24j |25 
liill this is truly classical. 

In what follows, we provide two proofs of a recursive scheme to obtain ex- 
plicit values of Rayleigh functions for a non-local BVP posed in UJ. Our main 
contribution is Theorems UJ [3j |4]for that non-local BVP. The proofs of the first 
two theorems are demonstrated directly using the properties of the eigenvalues 
of the non-local BVP without using the trace formulas unlike the way Goodwin 
proved for the regular BVPs 124]. The proof of Theorem |4] uses the generating 



functions as Radoux 120J and Liron [21] did for different BVPs (see also Ismail 



andMuldoon 115J). 



2. A Non-local Boundary Value Problem 

In this section we describe our eigenvalue problem. We first recall Corollary 
6 from the article LU : 
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Corollary 1. The eigenfunctions of the integral operatorX with the kernel K(x, y) = 
-\x- y\/2 for the unit interval Q = (0, 1) satisfy the following Laplacian eigen- 
value problem: 

-0" = A0, x£(0,l); 
0(O) + 0(1) = -0'(O) = 0'(1), (3) 
which can be solved explicitly as follows. 

• Xq « -5.756915 zs z7ie smallest (and the only negative) eigenvalue and is 
the solution of the following secular equation: 



coth 



(4) 



77ze corresponding eigenfunction is: 



<po(x) = C cosh \/-Ao 



where Co = \/2 



1 + 



sinh y'-Ao 
\/-Ao 



-1/2 



« 0.7812598 is a normalization constant 

\/-ho j 

tohave\\(p \\ L 2 in) = 1. 

• A2m-i = (2m - l) 2 7r 2 , m = 1, 2, . . ., and the corresponding eigenfunction is: 

<p2m-i(x) = \/2cos(2ra- \)nx. 
These are canonical cosines with odd modes. 

• X 2m , m = 1,2,..., is the solution of the secular equation: 

V &2m V ^2m 



cot- 

2 2 

and the corresponding eigenfunction is: 



(5) 



(f>2m(x) = C 2TO COSv / A 2 m|^-^j» 



where C 2m 



1 + 



sin \J A 2 n 
V ^2m 



-1/2 



zs a normalization constant. 
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Remark 1. We refer the reader to 1 1] for the motivation of considering such an 
integral operator X, the description of the higher dimensional versions, and 
a variety of applications. Here, however, we would like to point out our new 
interpretation of the above eigenvalue problem that was not explicitly stated 
in Q1 . The above problem turns out to be equivalent to the following problem 
defined for the whole real axis and then restricting the solutions to the unit 
interval fi. 

„ J Ai/a for x £ D.; 
lo forx£lR\n, 



with the continuity conditions at the boundary points: iff(0-) = y/(0+), if/'(Q—) = 
if/'{0+), y/{l-) = y/{l+), yr'il-) = Then, 0(x) in CorollarydJis IqW^W. 

Remark 2. The three cases of the eigenvalues in Corollary[T] i.e., A ; {X 2m -\}; 
and {Xzm} can also be derived from a single equation: 

( e <xl2 + e -al2 

-l = o, 



r a/2, -al2\ e +e a 
1 ' \ e al2_ e -al2 2 ) 



where A = -a 2 , and a £ C. Searching zeros of the first factor for a e iU leads to 
A2OT-1 - (2m- l) 2 n 2 whereas doing so in the second factor for a £ R leads to d4j 
and for a £ ilR leads to 10. 

Remark 3. Both Radoux ji^l and Liron |2ll] dealt with the secular equation 
tan f5 = f5. They explicitly mention that this equation came from the one-dimensional 
Laplacian eigenvalue problem by setting A = - a 2 , a = i/3, /5 £ R with the follow- 
ing Robin boundary condition: 

0(0) = 0, 4/(1) = 0(1). 

Note that 0'(O) = 0(0), 0(1) = lead to tan/3 = -j8. 

On the other hand, Radoux also dealt with the other secular equation cot/3 - 
f) whereas Liron treated the case involving cot/3 = -/3. Neither of them ex- 
plained why they wanted to treat these secular equations and neither of them 
explicitly listed the corresponding boundary condition unlike the case of tan/3 = 
p. In fact, simple computations similar to those in 10 Sec. 4.3] suggest that 
cot/3 = disassociated with the Robin boundary conditions (0'(O),0'(1)) = (0(0), 0) 
or (0, -0(1)) and cot/3 = -/3 is associated with (0'(O), 0'(1)) = (0, 0(1)) or (-0(0), 0). 
But one also needs to consider the hyperbolic versions, i.e., cotha = a, in or- 
der to fully solve the eigenvalue problems with the Robin boundary conditions 
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(0'(O),0'(1)) = (0,0(1)) or (-0(0), 0). Note that these Robin boundary condi- 
tions are all decoupled, i.e., local. To the best of our knowledge, [ 1] is the first to 
explicitly describe the unusual non-local boundary condition lf3j. 



Remark 4. One can exploit the well-known trace formula 1 24| |25| |26[ 

00 1 f 1 

E T? = / K p {x,x)dx, (6) 

n=0 A„ JO 

where K p {x, y) denotes the pth iterated kernel of K{x, y), to determine the first 
few expressions for the Rayleigh function at hand. Indeed, one obtains at once 



n =0 A n JO 



_ K{x,x)dx = 0, 

n =0 A n Jo 

and 

~ 1 f 1 , x x 2 \ , 1 

> -rr= I iC 2 (x,x)dx = / + — dx = — . 

tott Jo Jo U 4 4 J 24 

However this task becomes tedious for p > 3, and we propose to obtain these 
power sums without recourse to iterated kernels, but by exploiting properties 
of the transcendental equations of which the eigenvalues are roots. 



3. Sum of the Reciprocals of the Eigenvalues of CorollaryQ] 

In light of RemarkHl we want to show the following directly. 

Theorem 2. Let {A„}^1 be the eigenvalues of the boundary problem in Corol- 
lary^ andletK{x,y) = - \x-y\l2. Then, they satisfy the following trace formula: 

00 1 f 1 

E — = / *T(x,x)dx = 0. 
w =o A« Jo 

Proof. Let us group the eigenvalues into the three groups as indicated in Corol- 
lary!] 

S 1 1 S 1 SI 

y — = — + y + y — . 7 

n =0 A n A m =i A 2m-1 TO =i A 2m 
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Now, the second term of the sum is: 



E 



m -\ ^2m-l 



E 
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= al 



ti (2m-l)¥ 
1 



7T 



Vm=l 



1 (2m)' 



1 3 S 1 
Z2 " 7 2^ ~~a 



n 



m 1 - 



1 3 7l l 
7^ "i "T 



1 

8' 



(8) 



where we used the famous Basel problem identity L^ =1 1/ m 2 = tt 2 /6 resolved 
byEuler (see also HUH HQ). 
As for the last term of 0, 



E 



iA 



2m 



1 00 1 

4 ^ n x^' 

m= 1 m 



(9) 



where x m := \Z%im/2 > is the mth zero of the following transcendental equa- 
tion; see ©: 

cotx= -x. (10) 

To proceed to compute (9) explicitly, let us analyze QO) more deeply. Following 
Radoux 1 20] , let us first consider the following function and its Maclaurin series 
expansion: 



(cotx + x) -sinx 



cos x + xsinx 

1 — + — 

2! 4! 



v 

X 2 ( 1 

1 + — 

2 



+ x- 



X + 

13! 4! 



3 5 
x x 

X + — 

3! 5! 

15! 6! J 



1 + jc* + -x D 

2 4! 6! 



+ (-D 



ifc-l 



2fc-l 
(2*0! ' 



.2/c 



(ID 



+ ■ 



Now, the function cos x+ xsinx can also be expanded into the following infinite 
product in a manner similar to what Euler 1 17] and Rayleigh [5] did (see also 
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mi nana us, 26i): 



cos x + xsinx : 



1+3 n 



a 



00 ' x 2 



(12) 



where a « 1.19967864 satisfies a = cotha. 

In other words, x = ±ia are the two (and only) pure imaginary roots of 
cosx + xsinx. This can be verified as follows. Let us seek for the pure imagi- 
nary zeros of cos x + xsinx by setting x = iy, y eU. Then, we have 

cos(iy) + iysin(iy) 
e i(iy) + e -Hiy) gi(iy) 

- + iy — 

e y- e -y 



cos x + xsinx 



-i(iy) 



2 

e^ + e^ 



y- 



2i 
= 0, 



2 ' 2 
which is equivalent to coshy - y sinhy = 0, i.e., 

y = cothy. 

The justification for the product formula Q2J follows considerations similar to 
those for example in Gil Chap. 1]; see also I15L From Q2J, we have 



(13) 



cos x + xsinx 



m=l 
1 + 



n i-T ^n 



J2. 00 



(I 00 1 

Ot m-\ %m 



x 2 + . 



X 



(14) 



Equating the corresponding coefficients of the x 2 terms of CCD and CU, we 
have 

m-l A m ^ 

Hence, inserting this to IE), in turn, Q together with © gives us 

x J 1 00 1 cx) j 



^0 A n Ao ^1 A2to-1 7^1 ^2//; 

1 1 1 ~ 1 

= + _ + 

Ao 8 4 TO=1 x m 

1 11(1 1 

lo 8 4 la 2 " 2 

1 1 

A 4a 2 
= 0, 
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since A = -4a 2 , which can be verified by identifying 01 with the equation ([T3| 

via a = \/-X /2. □ 



4. Sums of Higher Powers of the Reciprocals of the Eigenvalues of CorollaryTT] 

Furthermore, we can establish the following identities: 

Theorem 3. Let {A„}^1 be the eigenvalues of the boundary value problem spec- 
ified in CorollaryUl Let K p {x, y) be the p th iterated kernel ofK{x, y) = -\x-y\l2. 
Then, we have 



OO l nl 

n=Q A-n JO 



1 I 

K p (x,x)dx = — 



(-D ? 

S2p + ~^F 



4 P_1 

+ \R 

2-{2p)V 



2p\ 



(15) 



where 



00 1 00 ( 4 
S 2p : = E = E 



777=1 X 



m-\ 



u 



2/;; 



andB 2p is the Bernoulli number, which is defined via the generating function: 

x _ y x n 
e x - 1 £o 

Moreover, Sz p satisfies the following recursion formula: 



"+ 1 {-l) n - e+1 {2{n- £ + !)-!) 
hi {2(n-£ + m 



S 2 e + 



t-lY 



a 



2f 



(-D" 

2(2 w)!' 



(16) 



Proof. The first equality in JT5T I connecting the sum of the powers of the eigen- 
values and the trace of the iterated kernel is the standard fact and its proof can 
be found in, e.g., |25[ Sec. 15]. Now, to prove the second equality, we have 



OO J 

E Tp 

n=0 A>n 



,p + E ,p 

A n m — i A" 



+ E 



L 



2771-1 mz 
| OO 



l A 2m 
1 



[4a 2 
(-D p 

4Pa 2 P ' 7T 2 p 
1 



+ 7T 2 p£ 1 (2m-l) 2 P ■ 4P^ 1X 2 P 



+ ■ 



1 n2n E 



1 5^ 1 

+ — E -^r- 

/I 77 ^ 2p 
111 

1 1 

+ 77t- S2 p 



2 2 p j ^ m 2 P 4P 



4^ S2 " + 



(-l)P 4^-1 ™ 1 

+ E 



4P 



2p a 2 P 
(-D p 



tt 2 p ^m 2 P 



S2p + 



+ 



4 P_1 



a 2 P j 2(2p) 
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where we used the following well-known formula first obtained by Euler (see, 
e.g., Illlll2lll3| to derive the last equality: 



{2n) 2 P 
l m 2 P 2{2p)\ 



00 i 

E 

m 



\B 2p \. 



Now, to prove the recursion formula i fTBl l. we follow Radoux |20] again. Tak- 
ing the logarithm of the product formula Q~2) followed by differentiation with 
respect to x, we have 



xcosx 



Zx 



cos x + xsinx 



V + E 



-2x 

v 2 



1 + ^2 

a 1 



which leads to 



-cosx = (cos x + x sin x) • < 
2 ^ 



00 1 1 

E 



" 1 H — 2 m -l m t 2" 



Expanding each term into the Maclaurin series or the geometric series, we have 



1 ~ (-l) w x 2w 



(2«)! 



E — — x 



(2fc)! 



E 



S 2 ^+2 I X 2£ 



2£+2 



Hence, comparing the coefficients of the x 2n term, we have: 



(-1)" 
2(2n)! 



n (_l)«-*-l(2«-2Jfc-l) 



E 



{2n-2k)\ 



(-1)' 



la 



2fc+2 



S2A: 



+2 



"+ 1 (-l)"-^ +1 (2(n-/ + l)-l) 



E 



(2(7»-/ + 1))! 



S 2 £ + 



i-iY 



a 



which is JTBT l. 



via setting f = fc+ 1, 

□ 



Let A p := L^l 7p • Here are me nrst f ew sums: 



1 1 

Ai = 0; A 2 = — ; A 3 = 

24 240 
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5. The Generating Function and Obtaining Recursive Formulas All at Once 



In this section, we show how to obtain the recursion formulas for the Ag 's at 
once and without recourse to the knowledge of Bernoulli numbers. The main 
result is the following theorem. 

Theorem 4. Let {X n }^L Q be the eigenvalues of the boundary problem in Corol- 
laryU\ and let K p {x, y) bethepth iterated kernel of K(x,y) = -\x-y\l2. Then, 



oo i r l 

n=0 A n JO 



K p (x,x)dx 



satisfies the recursion formula: 



n-l i 2 

4A M+1 + £(-l) fc 



fc=i 



{2k)\ {2k-\)\ 



A 



n-k+l 



{-\) n+l n 
(2w + 1)! ' 



n= 1,2,. 



with Ai = 0. 

Proof. From the statement of Corollary [T] and (HU , it is clear that 



(cosx + xsinx) -cosx = 



i+-2 n 



a 



where A = -4a 2 as defined above and where we set x^ = \f\kJ2, for k = 1, 2, . . .. 
One can again justify this product formula as in Knopp 1 18, Chap. 1] orany stan- 
dard Complex Analysis textbook which treats the Weierstrass Factor Theorem. 
In terms of the eigenvalues one has, after some trigonometric substitutions, 



1 + cosx x . 

+ — sinx ; 

2 4 



i-f n 



A, 



(17) 



Expanding the LHS into a Maclaurin series and equating lead to 



l+I(-l) fc 



fc=l 



J2k 



[2(2k)\ 4(2fc-l)! y 



£_L X 2 + £ _L_ 
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Witha fc :=(-l) fc (^ 



j denoting the coefficients of the Maclaurin ex- 



( 2(2k)\ 4(2fc-l)! 

pansion, one can recourse to Speigel's formulas Il9ll24ll 



OO 1 

k=0 A k 



- 



k=0 % 

OO J 

oo y 

E TT ~~ 



-ai 



■2a? 



3aia2 -3a3 - a. 



a\-Aa\a2 + 2a\ + Aaia2, -4a 4 



to obtain, as above, 



OO 1 

k=0 A k 

S 1 1 

to 4 24 
S 1 1 

^3 = E -» = — 

fc=o 4 240 
S 1 41 

^=Ett 



fc=0 A 4 fc 40320' 

One can generate a recursion formula for the At sequence employing what 
Ismail and Muldoon Qji call, properly, the "Euler-Rayleigh" technique. The 
logarithmic derivative of the entire function f(z) = 1+c ° sz + 1 sinz appearing in 
JT71 I gives, 

sinz zcosz 
"~ + ^~ _ 2z 2 ^ z 
1 + cosz zsinz A - z 2 A k ~ z 2 

2 4 

Or, substituting Ao = -4a 2 and A# = Ax 2 k , and after some manipulation, 



sin2z zcos2z 
+ 



l + cos2z zsin2z a 2 + z 2 £^ x 2 - z 2 
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£- r=:-zG(z). (18) 



The function 



Git) 



1 00 1 



« 2 +f 2 £i4- f2 



is known as the generating function of Ag. That is, one can obtain the needed 
recursion formula for this sequence from consideration of this function. To 
simplify notation, we let Mp := 4 £+1 Af + ±. It is then clear that 

(-1)' ~ 1 
Moreover, a straightforward calculation leads to 



Y,M e P e = Git). 



By dm, one then obtains 

sin2£ cos2£ fl + cos2f t 



At 



Expanding into power series leads to 



00 A n n I 00 1 n2k-\ 



n=l (2/1 + 1)! v fc=1 

From which one obtains M = 0, and 



j2Jfc— 2 



] t 2k \\YM e t 2e . 
[i2k)\ (2fc-l)!j ILtj, J 



/ 9 2fc-l 



)2fc-2 



(2jfc)! (2fc-l)! 
In terms of the Ag's one has, Ai = 0, as before, 



Me 



{-l) n+1 4 n n 
(2n + l)! 



4^ W+1 +E(-D 4 



>2fc-l 



>2fc-2 



I (2/c)! (2fc-l)!j (2ti + 1)I 

which is the same as the desired statement of the theorem. 



□ 



Remark 5. We note that the recursion generates the following valuesA 2 = 1/24, 
A 3 = -1/240, A 4 = 41/40320, A 5 = -107/725760, etc., corresponding to what we 
obtained differently in SectionlH 
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Remark 6. As in 1151, one can exploit the formulas generated for the A/s to 
obtain 

- I^m-ir 1 '^-" < A < -A- l ^ 2m] (19) 

and 

Mml Mm+\ < Ao < A 2 m-\l Mm- (20) 

for m = 1,2,3,.... These inequalities provide strict improvable bounds for the 
unique negative root of the transcendental equation 0J and another way of 
obtaining it. 

6. Higher Dimensional Considerations 

One of the motivations that led to the non-local BVP considered in jH is that 
one is able to read the spectral data (eigenvalues, eigenfunctions) by discretiz- 
ing then computing integrals involving the kernel K{x,y) over a domain D. c U d 
without imposing conditions on dQ. For the two dimensional case, K{x, y) takes 
the form of a logarithmic kernel 

K{x,y) = -— ln|x-y|. (21) 
in 



Troutman |28|] gave an analytical proof for the existence of at most one nega- 
tive eigenvalue and gave an upper bound estimate for it in terms of the area 
and transfinite diameter of H. (The transfinite diameter is a measure of the 



compactness of a domain; see 128|] for the definition.) In |29] , Kac offers a prob 



abilistic proof of this fact (see also 130J and the generalization in 131|]). With A p 
denoting the power sum in © and the iterated integrals computed numerically 
COD and (EUJ provide a practical and improvable means of computing this neg- 
ative eigenvalue for a specific domain. When the transfinite diameter of D, is 
less than or equal to one, this negative eigenvalue disappears. This is the case 
of the unit disk. In (H, it was found that the eigenvalues of the nonlocal BVP 
associated with the kernel lED are of two types, n , with multiplicity 3, and 
im-\ n ™ m multiplicity 2, for m = 2,3, . . ., and n = 1,2, . . .. Based on the values 
of a 2 ({'v), one can generate for the first few power sums. While T^ =i 1/A^ is 
easily seen to diverge, 

^1 °° 3 I In 2 3\ 

V — = 3a 4 (0)+2y a 4 (v) = — + . 

ti^\ v=i 32 8|6 2j 
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Similarly 

oo y 00 

£ —j = 3<T 2^ (0) + 2 ^ cr 2£ (v) 

fc=l \ v=l 

can be carried out explicitly for £ = 3,4,..., but there may not be an obvious 
recursion scheme. 
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